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We propose a dynamic model for a system consisting of self-propelled agents in which the influence 
of an agent on another agent is weighted by geographical distance. A parameter a is introduced to 
adjust the influence: the smaller value of a means that the closer neighbors have stronger influence 
on the moving direction. We find that there exists an optimal value of a, leading to the highest 
degree of direction consensus. The value of optimal a increases as the system size increases, while 
it decreases as the absolute velocity, the sensing radius and the noise amplitude increase. 
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I. INTRODUCTION 


The collective motion is a ubiquitous phenomenon in 
nature .examples of which include traffic jams 1], bird 
flocks [Ml, fish schools [f|, insects swarms [3, 0, bac¬ 
teria colonies [§, ftl , pedestrian flows EMI and active 
granular media [lj]. In recent years, a variety of efforts 
have been devoted to modeling the dynamic properties 
of swarms fl5l - f27l |. In 1995, Vicsek et al. proposed a 
particularly simple but rich model (28}. In the Vicsek 
model (VM), some self-propelled agents move with the 
same absolute velocity in a square-shaped cell with the 
periodic boundary conditions. At each time step, every 
agent updates its direction according to the average di¬ 
rection of the agents’ motion in its neighborhood. The 
neighborhood of an agent i is composed by agent i itself 
and those agents who fall in a circle of sensing radius 
that centered at the current position of i. It has been 
demonstrated that all agents will converge to the same 
direction on a macroscopic scale when the density of the 
system is high and the noise is small enough [291]. 


The VM and its variations have attracted much at¬ 
tention in the past decade j30l - [39| . Gregoire and Chate 
found that the onset of collective motion in the VM as 
well as in related models with and without cohesion is 
always discontinuous j3(3]. Hucpe and Aldana studied 
intermittency and clustering in the VM [3lj]- Yang et 
al. considered a power-law distribution of sensing ra¬ 
dius which can enhance the convergence efficiency [32j . 
Li et al. proposed an adaptive velocity model in which 
each agent not only adjusts its moving direction but also 
adjusts its speed according to the deg ree of direction con¬ 


sensus among its local neighbors [33| . Tian et al. discov¬ 


ered that there exists an optimal view angle, leading to 
the fastest direction consensus Q. Gao et al. proposed 


a restricted angle model that significantly improves the 
collective motion of self-propelled agents [33] . Schubring 
and Ohmann proposed a density-independent modifica¬ 
tion of the VM in which an agent interacts with neighbors 
defined by Delaunay triangulation j3l| . Peruani and Bar 
demonstrated that the clustering statistics and the corre¬ 
sponding phase transition to non-equilibrium clustering 
found in many experiments and simulation studies with 
self-propelled particles with alignment can be obtained 
by a simple kinetic model 1-1 fl . Buscarino et al. found 
that the direction consensus can be improved by long- 
range interactions [M EH ■ 

In the original VM, all agents in agent i’s neighbor¬ 
hood have the same influence on the moving direction 
of agent i. However, due to the existence of diversity in 
human society and animal world, influences of different 
agents usually are not the same. It has been shown that 
the heterogeneous influence plays an important role in 
various dynamics, such as the formation of public opin¬ 
ion S3 S3 and the evolution of cooperation S3 • In this 
paper, we propose a weighted VM, where the influence of 
neighbor j on agent i is determined by the geographical 
distance between the two agents. We set the weight of 
neighbors to be an exponential function with a tunable 
parameter a. Interestingly, we find that there exists an 
optimal value of a, leading to the highest degree of di¬ 
rection consensus. The effects of the moving speed, the 
sensing radius, the system size and the noise amplitude 
on direction consensus of the system are also studied. 

The paper is organized as follows. In Sec. El we intro¬ 
duce the distance-based influence model. Simulation and 
discussion are given in Sec. IIIII The paper is concluded 
in Sec. IIVI 


II. THE DISTANCE-BASED INFLUENCE 
MODEL 
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We consider N agents moving in the two-dimensional 
plane without periodic boundary conditions [33l libj . Ini¬ 
tially agents are randomly distributed on a region of L x L 
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rectangle with random directions. Note that this rectan¬ 
gle does not represent the boundary for motion, but only 
restricts the initial distribution of positions of agents. 
Each agent has the same absolute velocity vo and sensing 
radius r. At time t, the position of a specific agent i is 
updated according to 

Xi(t + 1) = Xj(t) +v 0 e l6 ^\ (1) 


Its direction is updated as 


e i0i(i+l) 


II E j6 r iuWjVeW II’ 


( 2 ) 


where A 6i € [— 77 , rj\ denotes the white noise, e 19 ^ is 
a unit directional vector, Ti(t) is the set of neighbors of 
agent i defined by the sensing radius r at time step t, and 
the weight Wj (t) denotes the influence of the neighbor j 
of agent i at time step t. We define Wj(t) as 

Wj(t) = e “- d « (t) , (3) 


where dij(t) = ||xj(i) — Xj(i)|| is the distance between 
agent i and j at time step t, and a is a tunable parameter. 
For a > 0(< 0), the farther (closer) neighbors have larger 
weight of influence. When a = 0, our model is the same 
as the standard VM, where each agent in the system has 
the same weight. 


III. SIMULATION AND DISCUSSION 

In all the following simulations, we set L = 10. 

We first consider the case in which the noise is zero 
(r) = 0). From the perspective of complex network the¬ 
ory, the topology of self-propelled agent system can be 
expressed as a temporal network [47| . At time t , each 
agent is represented by a node and an edge between agent 
i and j is established if the distance between them is 
shorter than the sensing radius r. A cluster is a subgraph 
in which any two nodes are connected to each other by 
paths running along edges of the network. In the case of 
zero noise, after a period of evolution, agents aggregate 
into different moving polar clusters. All agents within 
a moving polar cluster move in the same direction, as 
shown in Fig. [0 

Following the previous studies 0,0 , we examine the 
cluster size distribution. The results are shown in Fig. [2] 
We see that the probability P(S) that a cluster with size 
S decays as a power law, following P(S) ~ S ~U The 
inset of Fig. [2] shows the dependence of /i on a when 
the absolute velocity vo = 1 and the sensing radius r = 
1.6. We find that the exponent n is minimum at a « 4. 
Since smaller 77 indicates a higher probability of large 
clusters, a ~ 4 is an optimal value for the system to have 
large clusters. The value of /r is found to be in the range 
0.8 < n < 1.5 in line with the previous experiments and 
simulations 00,0- 



FIG. 1: Snapshots of locations and velocities (a) in the initial 
configuration, and (b) at the time step t = 100. The arrows 
show the direction of motion of the agents. The parameters 
are set as L = 10, N = 50, r = 2, vo = 0.2, and a = 0. 
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FIG. 2: (Color online) The distribution of cluster sizes P(S) 
for different values of a. It is found that P(S) ~ The 

inset shows that the exponent /r as a function of a. The abso¬ 
lute velocity vo = 1, the sensing radius r = 1.6 and the system 
size N = 1000. Each data point results from an average over 
10 4 different realizations. 


Following the previous studies [32;, 0,0, we measure 
the degree of direction consensus by the relative size Si 
of the largest cluster, which is defined as the ratio of the 
number of agents within the largest cluster to the total 
number of agents. Note that 0 < Si < 1. A larger value 
of Si indicates a higher degree of direction consensus. 

Figure [3] shows that si as a function of the absolute 
velocity vo for different values of a. From Fig. [3l we can 
see that for any given value of a, si decreases monoton- 
ically as vq increases. Figure [I] plots si as a function of 
the sensing radius r for different values of a. In Fig. [4J 
one can observe that for any given value of a, Si increases 
with r. Figure [5] depicts Si as a function of the system 
size N for different values of a. From Fig.[5j we find that 
for any given value of a, Si increases with N. 

Figure O shows the dependence of si on a for differ- 
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FIG. 3: (Color online) The relative size si of the largest clus¬ 
ter in the steady state, as a function of the absolute velocity 
vq for different values of a. The sensing radius r = 1.6 and 
the system size N = 1000. Each data point results from an 
average over 1000 different realizations. 



FIG. 4: (Color online) The relative size si of the largest clus¬ 
ter in the steady state, as a function of the sensing radius r 
for different values of a. The absolute velocity vo = 0.6 and 
the system size N = 1000. Each data point results from an 
average over 1000 different realizations. 


ent values of the absolute velocity vg- In Fig. EH one can 
observe a resonance-like behavior. For a given value of 
vo, there exists an optimal value of a, hereafter denoted 
by cx op t, leading to the maximum si- It is interesting to 
note that, for the absolute velocity vq = 1 and the sens¬ 
ing radius r = 1.6, a op t ~ 4 not only corresponds to the 
maximum si, but also results in the minimum /i (see the 
inset of Fig. n. Intuitively, at this point, a minimum n 
corresponds to a distribution P(S) that decreases most 
slowly as S increases, indicating a larger probability that 
a large cluster emerges. This infers that the minimum 



N 

FIG. 5: (Color online) The relative size si of the largest clus¬ 
ter in the steady state, as a function of the system size N 
for different values of a. The absolute velocity vo = 1 and 
the sensing radius r = 1.6. Each data point results from an 
average over 1000 different realizations. 



FIG. 6: (Color online) The relative size si of the largest clus¬ 
ter in the steady state, as a function of a for different values 
of the absolute velocity vo- The inset shows the optimal value 
of a, a op t, as a function of vo . The sensing radius r = 1.6 and 
the system size N = 1000. Each data point results from an 
average over 1000 different realizations. 


of /i and the maximum of Si occur at the same point of 
a opt ~ 4. From the inset of Fig. (6j we find that a opt de¬ 
creases from 10 to 4 as vo increases from 0.3 to 1. Since 
larger a indicates more influence from distant nodes, this 
means that as vg increases, direction consensus is bet¬ 
ter preserved by strengthening the influence from closer 
agents. 

Figure [7] shows the dependence of si on a for different 
values of the sensing radius r. From Fig. [71 we find that 
for a given value of r, there exists an optimal value of a 
leading to the maximum si- The inset of Fig. 0 displays 











4 



FIG. 7: (Color online) The relative size si of the largest clus¬ 
ter in the steady state, as a function of a for different values 
of the sensing radius r. The inset shows the optimal value of 
a, a op t, as a function of r. The absolute velocity vo = 0.6 
and the system size N = 1000. Each data point results from 
an average over 1000 different realizations. 



FIG. 8: (Color online) The relative size si of the largest clus¬ 
ter in the steady state, as a function of a for different values 
of the system size N. The inset shows the optimal value of 
a, a op t, as a function of N. The absolute velocity vo = 1 and 
the sensing radius r = 1.6. Each data point results from an 
average over 1000 different realizations. 


that a op t decreases from 8 to 3 as r increases from 1.2 to 
2.4. Figure H] shows the dependence of si on a for differ¬ 
ent values of the system size N. From Fig. [8j one can see 
that for a given value of N, there exists an optimal value 
of a resulting in the maximum si. The inset of Fig. [8] 
displays that a op t increases from 2 to 5 as N increases 
from 200 to 10000. 

The resonance-like behavior of the dependence of Si 
on a can be understood as follows. For very small values 
of a, an agent is mostly influenced by its closest neigh¬ 



FIG. 9: (Color online) The relative size si of the largest clus¬ 
ter as a function of the time step t for different values of a. 
The absolute velocity vo = 0.6, the sensing radius r — 1.6, the 
system size N = 1000 and the noise amplitude rj = 0.6. Each 
curve results from an average over 1000 different realizations. 


bors and the active interaction radius is much shorter 
than the sensing radius r. As indicated by the a = — 2 
case in Fig. [2 si increases very slowly with r. Previous 
studies have shown that the direction consensus can be 
improved by adding long-range interactions |U HoJ . 
Thus it is beneficial if distant neighbors are given a strong 
influence. However, when a is too large that an agent is 
mostly influenced by the farthest neighbor, the interac¬ 
tion between other neighbors would be so weak that is 
not able to maintain the coherent motion between them, 
which in turn also disrupts direction consensus. A sim¬ 
ilar effect has been observed in percolation issues [5l| . 
Therefore, there must be an intermediate value of a that 
optimizes direction consensus. 

Next, we study the case in which the noise amplitude 
r) is nonzero. Figure [9] shows that s± as a function of the 
time step t for different values of a when rj = 0.6. From 
Fig. 0 we see that si decreases to 0 as time evolves, 
indicating that all the agents will disperse away without 
any apparent cluster in the noisy environment. From 
Fig. M one can also observe that the value of a affects 
the process of dispersion. 

To quantify the speed of dispersion in the noisy case, 
we study the transient time r, which is defined as the 
time when si first below a certain value e. We have 
checked that qualitative results are invariant when e is 
small enough. In this paper we take e = 0.02. Figure [TO] 
shows the transient time r as a function of a for different 
values of the noise amplitude rj. One can see that r can 
be maximized by an optimal value of a. This means the 
process of dispersion can be best slowed down by fine 
tuning the value of a. Moreover, the inset of Fig. [TO] 
displays that the optimal value of a, a op t decreases from 
3 to 1 as the noise amplitude rj increases from 0.2 to 0.7. 

In the above studies, the degree of direction consensus 
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FIG. 10: (Color online) The transient time r as a function 
of a for different values of the noise amplitude 77 . The inset 
displays the optimal value of a, a op t, as a function of 77 . The 
absolute velocity Vo = 0 . 6 , the sensing radius r = 1.6 and 
the system size N = 1000. Each data point results from an 
average over 1000 different realizations. 



FIG. 11: (Color online) The order parameter V a as a function 
of a for different values of the noise amplitude 77 . The absolute 
velocity uo = 0 . 6 , the sensing radius r = 1.6 and the system 
size N = 1000. Here V a is calculated at the time step t = 500. 
Each data point results from an average over 1000 different 
realizations. 


of a system is quantified by the relative size of the largest 
cluster. We have checked that the qualitative results hold 
unchanged if we use an order parameter to measure the 
degree of direction consensus, as defined by [28} 

1 N 

V a = N H’ Q<Va<l. (4) 

i= 1 

A larger value of V a indicates a better consensus. Fig¬ 


ure EH shows the dependence of V 0 on a for different 
values of the noise amplitude 17 . From Fig. 111! one can 
see that for a given value of 77 , V„ is maximized by an 
optimal value of a. Besides, we find that V a decreases as 
77 increases when a is fixed. 


IV. CONCLUSION 

In conclusion, we have introduced a weight based on 
the geographical distance into the original Vicsek model. 
A tunable parameter a is used to govern the weight. The 
direction of each agent is updated by the weighted av¬ 
erage directions of its neighbors. In the noiseless case, 
self-propelled agents will aggregate into different clusters 
and the size distribution of clusters exhibits a power-law 
form. We find that direction consensus is enhanced when 
the absolute velocity is small, the sensing radius or the 
system size is large. Interestingly, there exists an opti¬ 
mal value of a that yields the highest level of direction 
consensus. In the noisy case, agents gradually disperse 
away and each agent moves lonely. This process of dis¬ 
persion can be best slowed down by optimizing a. The 
value of optimal a depends on the absolute velocity, the 
sensing radius, the system size and the noise amplitude. 
For all kinds of cases, the optimal value of a is always 
positive, indicating that a suitably strong influence of dis¬ 
tant neighbors can best enhance the direction conseusus 
of the system. 

Note that recently Gao et al. proposed another 
weighted Vicsek model [46}, where the weight of each 
agent is determined by its number of neighbors. They 
found that the system becomes easier to achieve direction 
consensus when agents with more neighbors are allocated 
greater weight. In the weighted model proposed by Gao 
et al ., agents are assumed to be able to obtain global in¬ 
formation, i.e., an agent not only knows the number of 
its neighbors, but also knows the number of neighbors’ 
neighbors. In our model, an agent only gets local infor¬ 
mation, i.e., the distance between it and its neighbor. 
All together Ref. [46[ and our work provide a deeper un¬ 
derstanding of the impact of heterogeneous influence on 
collective motion. 
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